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Siegel disks are domains around fixed points of holomorphic maps in which the maps are locally 
linearizable (i.e., become a rotation under an appropriate change of coordinates which is analytic in 
a neighborhood of the origin). The dynamical behavior of the iterates of the map on the boundary of 
the Siegel disk exhibits strong scaling properties which have been intensively studied in the physical 
and mathematical literature. 

In the cases we study, the boundary of the Siegel disk is a Jordan curve containing a critical point of 
the map (we consider critical maps of different orders), and there exists a natural parameterization 
which transforms the dynamics on the boundary into a rotation. We compute numerically this 
parameterization and use methods of harmonic analysis to compute the global Holder regularity of 
the parameterization for different maps and rotation numbers. 

We obtain that the regularity of the boundaries and the scaling exponents are universal numbers 
in the sense of renormalization theory (i.e., they do not depend on the map when the map ranges in 
an open set), and only depend on the order of the critical point of the map in the boundary of the 
Siegel disk and the tail of the continued function expansion of the rotation number. We also discuss 
some possible relations between the regularity of the parameterization of the boundaries and the 
corresponding scaling exponents. 

PACS numbers: 05.45.Df, 02.30. Nw, 05.10.-a, 05.10.Cc 



According to a celebrated theorem by Siegel, 
under certain arithmetic conditions, the dynami- 
cal behavior of the iterates of a holomorphic map 
around a fixed point of the map is very simple 
the iterates of the map fill densely analytic topo- 
logical circles around the critical point. In the 
domain around the critical behavior where the it- 
erates exhibit such behavior called Siegel disks, 

there exists a complex analytic change of vari- 
ables that makes the map locally a multiplica- 
tion by a complex number of modulus 1. On the 
boundary of a Siegel disk, however, the dynami- 
cal behavior of the iterates of the map is dramati- 
cally different for example, the boundary is not 
a smooth curve. The iterates on the boundary 
of the Siegel disk exhibit scaling properties that 
have motivated the development of a renormaliza- 
tion group description. The dynamically natural 
parameterization of the boundary of the Siegel 
disk has low regularity. We compute accurately 
the natural parameterization of the boundary and 
apply methods from Harmonic Analysis to com- 
pute the Holder exponents of the parameteriza- 
tions of the boundaries for diffterent maps with 
different rotation numbers and with different or- 
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ders of criticality of their critical points. 

I. INTRODUCTION 

Siegel disks - the domains around fixed points of holo- 
morphic maps in which the map is locally linearizable 
(defined in more detail in Section [ll]) - are among the 
main objects of interest in the dynamics of holomorphic 
maps. Their boundaries have surprising geometric prop- 
erties which have attracted the attention of both math- 
ematicians and physicists. Notably, it was discovered 
in 1,2] that in some cases there were scaling relations 
for the orbit, which suggested that the boundary was a 
fractal object. Since then this phenomenon has been a 
subject of extensive numerical and mathematical studies 
[Slll|51ISl|71[il[Tni[IIJ[Tl]. 

In this paper, we report some direct numerical calcu- 
lations of the Holder regularity of these boundaries for 
different rotation numbers of bounded type and for dif- 
ferent maps. 

The main conclusion of the numerical calculations in 
this paper, is that, for the cases we consider, the bound- 
aries of the Siegel disks are C curves for some k > 0, 
and we can compute numerically the value of k. Even if 
we - obviously ~ consider only a finite number of cases, 
we expect that the results are significative for the Siegel 
disks of polynomials with rotation numbers which have 
an eventually periodic continued fraction. 

The values of the Holder regularity k are, up to the 
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error of our computation, universal in the sense of renor- 
malization group analysis, namely that they are indepen- 
dent of the map in a small neighborhood in the space of 
maps. We also performed computations for maps whose 
rotation numbers have the same "tail" of their continued 
fraction expansion, and found that our numerical results 
depend only on the tail. 

Our computation of the Holder regularity are based 
on the method introduced in [T3], which is a numerical 
implementation of several constructions in Littlewood- 
Paley theory. This method was also used in [TH |TH1 |T5J 
HZ]. 

For the case of the golden mean rotation number, the 
fact that the boundaries of Siegel disks are Holder was 
proved in assuming the existence of the fixed point of 
the renormalization operator conjectured in [2\. The ex- 
istence of a fixed point of a slightly different (and presum- 
ably equivalent) renormalization operator was proved 
in It seems that a similar argument will work for 
other rotation numbers with periodic continued fraction 
expansion provided that one has a fixed point of the 
appropriate renormalization operator. These arguments 
provide bounds to the Holder regularity n, based on prop- 
erties of the fixed point of the renormalization operator. 

Our computations rely on several rigorous mathemat- 
ical results in complex dynamics. Notably, we will use 
that that for bounded type rotatation numbers and poly- 
nomial maps, the boundary of the Siegel disk is a Jordan 
curve, and contains a critical point [TH]. See also [20] 
for other results in this direction. We recall that a num- 
ber is bounded type means that the entries in the contin- 
ued fraction of this number are bounded. Equivalently, 
a number a is of constant type if and only if for every 
natural n and integer m we have \an — m\^^ < vn~^ [2T] . 

It is also known in the mathematical literature that 
for non-Diophantinc rotation numbers (a case we do not 
consider here and which indeed seems out of reach of 
present numerical experiments), it is possible to make 
the circle smooth [251 HSl 121] or, on the contrary, not a 
Jordan curve [25| [25] . 

The plan of this paper is the following. In Section [H] 
we give some background on Siegel disks and explain how 
to parameterize their boundaries. Section |III| is devoted 
to the numerical methods used and the results on the 
regularity of the boundaries of the Siegel disks and the 
scaling properties of the iterates. In Section [TV] we con- 
sider some connections with geometric characteristics of 
the Siegel disk (in particular, its area), in Section [v] we 
derive an upper bound on the regularity, and in the final 
Section IVlj we recapitulate our results. 



II. SIEGEL DISKS AND THEIR BOUNDARIES 
A. Some results from complex dynamics 

In this section we summarize some facts from complex 
dynamics, referring the reader to [27] [25] for more details. 



We consider holomorphic maps of C that have a fixed 
point, and study their behavior around this point. With- 
out loss of generality, in this section we assume that the 
fixed point is the origin, so that the maps have the form 



f{z) = az + 0{z^) 



(1) 



(For numerical purposes, we may find more efficient to 
use another normalization.) We are interested in the case 



that \a\ = 1, i.e., a 



where a € [0, 1) is called the 



rotation number of /. In our case, we take / to be a 
polynomial, so that the domain of definition of the map 
/ is not an issue. 

The stability properties of the fixed point depend cru- 
cially on the arithmetic properties of a. The celebrated 
Siegel's Theorem [2211301131] guarantees that, if ct satis- 
fies some arithmetic properties (Diophantine conditions), 
then there is a unique analytic mapping h (called "con- 
jugacy" ) from an open disk of radius r (called the Siegel 
radius) around the origin, i?(0, r), to C in such a way 
that h{0) = 0, h'{0) = 1, and 



/ o h(z) = h(az) . 



(2) 



We note that the Siegel radius is a geometric property 
of the Siegel disk. It is shown in [32] that h can be 
characterized as the conformal mapping from B{0,rs) to 
the Siegel disk mapping to and having derivative 1. 
Later in Section [IV A| we will show how the Siegel disk 
can be computed effectively in the cases we consider. 

We refer to f331 [31] for some mathematical devel- 
opments on improving the arithmetic conditions of the 
Siegel theorem. In this paper we only consider rotation 
numbers that satisfy the strongest possible Diophantine 
properties. Namely, we assume that a is of bounded type 
and, in particular, perform our computations for num- 
bers with eventually periodic continued fraction expan- 
sions (see Section II C for definitions). In this case, there 
is an elementary proof of Siegel's Theorem [35". 

Let rg stand for the radius of the largest disk for which 
the map h exists. The image under h of the open disk 
B{0,rs) is called the Siegel disk, V, of the map /. For 
r < Ts, the image of each circle {w € C : \w\ — r} 
under h is an analytic circle. The boundary, dT>, of the 
Siegel disk, however, is not a smooth curve for the cases 
considered here. 

The paper ^ contains numerical observations that 
suggest that the dynamics of the map / on dV satisfies 
some scaling properties. These scaling properties were 
explained in certain cases by renormalization group anal- 
ysis [21 HI HH [m [36]. These scaling properties suggest 
that the boundaries of Siegel disks can be very interesting 
fractal objects. 

Clearly, the Siegel disk cannot contain critical points 
of the map /. It was conjectured in [Tj that the boundary 
of the Siegel disk contains a critical point. The existence 
of critical points on the boundary depends on the arith- 
metic properties of a and it may be false [25] , but it is 
true under the condition that the rotation number is of 
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FIG. 1: Critical points c and Siegel disks of the map Q 
for cr = (1°°) = 1(^5 - 1) (* and a thin hne), a = (5°°) = 
1(^29-5) (+ and a thick line), and a = (51°°) = ^(9-%/5) 
(x and a dotted line). 



bounded type [TB], which is the case we consider in this 
paper. See also [37] . 

In Figure [T] we show the Siegel disks of the map 



'z + zJ 



(3) 



for different rotation numbers a (for the notations for 
a see Section II C). In all cases the only critical point, 

c : 



1 ^27riCT 
2 



e— , is simple: /'(c) = 0, /"(c) = 2^Q. 



In this paper we study maps of the form 
1 



[5„,+i(z)-(l-/3)5™(z)] , (4) 



where m e N, /3 is a complex parameter, and the function 
Qm : C ^ C is defined as 



1 



-[l-(l-z)™+i] 



Let / be a map of the form ([ij, and c be its criti- 
cal point that belongs to the boundary of the Siegel disk 
of this map (we will only consider cases where dT) con- 
tains one critical point). Let d be the multiplicity of the 
critical point c, i.e., f'^^\c) = for /c = 1, 2, . . . , d, and 
^ 0. We will call d the order of the critical 

point. 

Noticing that, for the map Q, 



fL.A^) = (1 - zy 



1 - 



/3 



we see that /m,a,/3(0) = 0, /^^^,^(0) = e^"""^, and, more 
importantly, \l (5 ^ 1 , the point z = 1 is a zero of /^^ ^ ^ 




FIG. 2: Siegel disks of the maps /d,(2°°>,i+3i (given by Q) 
with critical point c = 1 of order d = 1 (outermost curve), 
d = 5, and d = 20 (innermost curve) - see text for details. 



of multiplicity to, while, for /3 = 1, the point z = 1 is 
a zero of fl^^p of multiplicity m + 1 . As long as the 
critical point z = (3 is outside the closure of the Siegel 
disk, the scaling properties of the iterates on dT> in the 
case of Diophantine cr are determined by the order d of 
the critical point z ^ \. Below by "critical point" we 
will mean the critical point that belongs to dT>. One of 
the goals of this paper is to study how the regularity and 
the scaling properties depend on the order of the critical 
point. 

In Figure |2] we show about 16 million iterates of the 
critical point c = 1 of the maps /d.(2°°),i+3i for order 
d = 1, 5, and 20, of the critical point c = 1 (the other 
critical point, [3 = l-l-3z, is not in the closure of the Siegel 
disks, so is irrelevant for the problem studied). Note that, 
especially for highly critical maps, the iterates approach 
the critical point very slowly because the modulus of the 
scaling exponent becomes close to 1 (see Table [ll|. 



B. Parameterization of the boundary of a Siegel 
disk 

In the cases considered here, the boundaries of Siegel 
disks cannot be written in polar coordinates as r = R{9) 
(because some rays 9 = const intersect dV more than 
once). In this section we explain how to parameterize 
dV, and define the functions whose regularity we study 
numerically. 

Once we know that a critical point c is in the boundary 
of the Siegel disk V (which in the cases we consider is 
guaranteed by the results of [H]), it is easy to obtain a 
parameterization of the boundary which semiconjugates 
the map / to a rotation. 

It is known from the mathematical theory that h — 
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which is univalent in the open disk B{0, rs) - can be ex- 
tended to the boundary of 5(0, rs) as a continuous func- 
tion thanks to the Osgood- Taylor-Caratheodory Theo- 
rem (see, e.g., [SS', Section 16.3] or [351 Section IX. 4]). 

A dynamically natural parameterization x of dV is ob- 
tained by setting 



x(t) = /.(rse2"(*+«)) , 



(5) 



where 6* is a constant to be chosen later. From it 
follows that 



/ o x{t) = x{t + (t) 



(6) 



with x(0) — c. Since we know that in our cases the 
critical point c is in dV, \h~^{c) | = rg, and we can choose 
9 so that 



x(0) = M^se^"^'')=c. 



(7) 



In summary, the function x : IT — > C defined by (|5| and 
([t]) is a parameterization of dV such that the dynamics 
of the map / on dV is a rotation by a on the circle T, 
and x(0) is the critical point. 

Iterating n times ([6]) for t — and using (|7]), we obtain 



X{na) = ric) 



(8) 



Since the rotation number tr of / is irrational, the num- 
bers na (taken mod 1) are dense on the circle T, and the 
iterates /"(c) of the critical point c are dense on dV as 
well. Hence, using equation ([8|, we can compute a large 
number of values of x by simply iterating /. 

Our main object of interest are the real, 5Rxi S'Hd imag- 
inary, 3xi parts of the function x- In Figure |3] we have 
shown the graphs of the real and imaginary parts of x pa- 
rameterizing the boundary of the Siegel disk of the map 
([3]) for two different rotation numbers. 



C. Continued fraction expansions and rational 
appr oximat ions 

In this section, we collect some of the notation on con- 
tinued fraction expansions. 

Let A = (ai, a2, . . . , ftp) be a finite sequence of p 
natural numbers aj € N; for brevity, we will usually 
omit the commas and write A ~ {aia2...ap). Let 
B — {bih2 ■ . - bq) be another sequence of natural num- 
bers bj e N, j = 1, 2, . . . , g, AB :— (0102 . . . 0^6162 ■ ■ - bq) 
stand for the concatenation of A and B, and A" stand for 
AA • • • A {n times). Let |A| = p denote the length of A. 

For a e N define the function Ta : (0, 1) (0, 1) by 

Ta{x) :— . Similarly, for A — (0102... «„), define 

the function JFa : (0,1) (0,1) as the composition 




FIG. 3; Graphs of the real (solid lines) and the imaginary 
(dashed lines) parts of the map x corresponding to dT> of the 
map ([3| for (J = (1°°) (thin lines) and a = (5°°) (thick hnes) 
- cf . Figure [l] 



Let (B) stand for the number whose continued fraction 
expansion (CFE) is given by the numbers in the se- 
quence B: 



(B) = (6162... fe,) := 



61 



'^a^x) 



the numbers hj are called the (partial) quotients of (B). 
A number a — (aia2 . . .) is of bounded type if all numbers 
Oi (z e N) are bounded above by some constant M. 

We are especially interested in studying numbers with 
CFEs of the form 

(AB°°) lim (AB") , 

n — >oo 

which are called eventually periodic (or preperiodic) . 
Since each number of this type is a root of a quadratic 
equation with integer coefficients (see [121 Theorem 176]), 
such numbers are also called quadratic irrationals. We 
will call A the head and B°° the tail, B the period, and 
|B| the length of the period of the CFE. 

If two quadratic irrationals have the same tail, they 
are said to be equivalent. One can prove that cr and p are 
equivalent if and only if cr = ('-P+ A) /{fip+ v), where the 
integers t. A, /i and v satisfy lv — Xfi — ±1 [40, Theorem 
175]. 



D. Scaling exponents 

Let / be a map of the form ([T]) with an eventually 
periodic rotation number a = (AB°°) with length of its 
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period q = |B|, and let c be the critical point of / on dV 
(and there are no other critical points on dT>). Let 



Qn 



(AB'= 



(9) 



where and Qm are natural numbers that have no 
common factors. Define the scaling exponent 



a 



lim 

rn — >oo 



/Q™(c)- 



(10) 



This exponent is a complex number that depends on the 
tail B of the CFE of a and on the order of the critical 
point c, but does not depend on the head A or on details 
about the map /. 

If the length |B| of the tail B of the CFE of the rotation 
number of the map is more than 1 , then B is determined 
only up to a cyclic permutation, and the argument of the 
complex number a is different for different choices. That 
is why we give our data only for |a|. 



III. HOLDER REGULARITY AND SCALING 
PROPERTIES OF THE BOUNDARIES OF 
SIEGEL DISKS ~ NUMERICAL METHODS AND 
RESULTS 

A. Some results from harmonic analysis 



Let K — n + k' , where n e {0, 
We say that a function <j) : T - 



,2,...}, and k' e (0,1). 
T has (global) Holder 



regularity k and write G C"^(T) if k = n + k' is the 
largest number for which exists and for some con- 
stant C > satisfies 

^(")(t)„0(")(s) <C It-sl"' Vi,sGT. 

We call attention to the fact that we do not allow k to be 
an integer since otherwise the definition of Holder so that 
the characterizations we discuss later must be modified. 
In our problem, k turns out to be non-integer, so that 
the characterizations we discuss apply. 

In the mathematical literature, there are many char- 
acterizations of the Holder regularity of functions. Some 
surveys that we have found useful are [HJ • 

In the paper |13j . we developed implementations of 
several criteria for determining Holder regularity numer- 
ically based on harmonic analysis, and assessed the reli- 
ability of these criteria. In this paper, we only use one of 
them, namely the method that we called the Continuous 
Littlewood-Paley (CLP) method which has been used in 
[Tl[T5J[Tni[T7]. The CLP method is based on the follow- 
ing theorem (tlj |42] : 

Theorem III.l. A function (p £ C"^(T) if and only if 
for some 77 > there exists a constant C > such that 
for all T > and rj > 



d_ 



(11) 



where A stands for the Laplacian: A0(t) = 4>"(t). 



Note that one of the consequences of Theorem |III.1| is 
that if the bounds (111 hold for some 77, they hold for any 
other r/. Even if from the mathematical point of view, all 
values of 77 would give the same result, it is convenient 
from the numerical point of view to use several to assess 
the reliability of the method. 



B. Algorithms used 



The 



algorithm we use is based on the fact that 
-rpp J c is a diagonal operator when acting on a 

Fourier representation of the function: if 



(j){t) = '(pk e 
fceN 



then 
dr 



(-27r|fc|)''e- 



'2Trr|fc| 2 „-27ri/ct 



The Fourier transform of the function (f) can be com- 
puted efficiently if we are given the values of (/)(i) on a 
dyadic grid, i.e., at the points t„i = 2^^^m, where M is 
some natural number and m ^ 0,1, . . . , 2^^ — 1. Unfortu- 
nately, the computation indicated in ([s]) gives the values 
of the function x = + * on the set {na}'^^Q of 
translations by the irrational number a. Therefore, we 
need to perform some interpolation to find the approxi- 
mate values of x on the dyadic grid, after which we Fast 
Fourier Transform (FFT) can be computed efficiently. 

Hence, the algorithm to assess the regularity is the 
following. 

1. Locate the critical point c (such that /'(c) = 0). 

2. Use equation ^ to obtain the values of the function 
X at the points {na}'^^Q for some large TV. 

3. Interpolate 3?x and 9% to find their approximate 



values on the dyadic grid {2 m} 



lm=0 ■ 



4. For a fixed value of rj, compute the value of the left 
hand side of (111 for several values of r by using 
FFT (for 



Sftx and separately for 
this for several values of rj. 



5x); do 



L-(T) 



5. Fit the decay predicted by (111 to find the regular- 
ity K. 

Let us estimate the cost in time and storage of the 
algorithm above keeping 2*^ values of the function. Of 
course, locating the critical point c is trivial. Iteration 
and interpolation require 0{2^) operations. Then, each 
of the calculations of (11) requires two FFT, which is 
0(2^1n(2^)) = 0{M2^). In the computers we used 
(with about 1 GB of memory) the limiting factor was 
the storage, but keeping several million iterates of / and 



FIG. 4: Plot of login I (^)fcl vs. log^j, |fc| for the map Q FIG. 5: Plot of log^, ] (Kx)fc| vs. logn, \k\ for the map Q 
with a = (1°°) - see text for details. with a = (5°°) - see text for details. 



computing 2^^ « 8x10^ Fourier coefficients of 5Rx and 5^ 
was quite feasible. (Note that a double precision array of 
2^^ double complex numbers takes 2^'' bytes — 128MB 
and one needs to have several copies.) 

The iterates f"'{c) were computed by using extended 
precision with GMP - an arbitrary-precision extension 
of C language |13]- This extra precision is very useful to 
avoid that the orbit scales. Note that the critical points 
are at the boundary of the domains of stability, so that 
they are moderately unstable. 

The extended precision is vitally important in com- 
puting the scaling exponents a. To obtain each value in 
Table [ll] we computed several billion iterates of the crit- 
ical point of the map. To reduce the numerical error, we 
used several hundred digits of precision. 



C. Visual explorations 

In Figures |4] and [s] we have plotted (with impulses) 
the modulus of the Fourier coefficients (3?x)fc versus |fc| 
on a log-log scale (for several million values of k) for the 
boundary of the Siegel disk corresponding to the map 
([3]) for rotation numbers a equal to (1°°) and (5°°), and 
order of the critical point d = 1 (the same cases as the 
ones in Figures [l] and |3]) . The self-similar structure of 
the boundary of the Siegel disk is especially clearly vis- 
ible in the "straightened-out" graph of the spectrum - 
in Figure [6]we plotted log^o \k (Jfix)k\ vs. log^p |fc| for the 
same spectra as in Figures |4] and |5] The width of each of 
the periodically repeating "windows" in the figure (i.e., 
the distance between two adjacent high peaks) is approx- 
imately equal to | logj^Q a\, where a is the corresponding 
rotation number. The periodicity in the Fourier series 
has been related to some renormalization group analysis 
in phase space [44j . 




FIG. 6: Plot of logi„ 
withCT = (1°°) (top), a 



^(5Rx)fcl vs. logj^„|fe| for the map (|3| 
= (5°°) (bottom) - see text for details. 



To illustrate the effect of the order of criticality on 
the Fourier spectrum of 3ffxj we showed in Figure ^the 
"straightened-out" graphs of the spectra, log;iQ |/c(3?x)fel 
versus logj^g |fc|, of the function SRx corresponding to the 
maps fd,{2'^),i+3i for orders d = 1,5,20 (the Siegel disks 
of these maps were shown in Figure [2]) . For all plots in 
this figure we used the same scale in vertical direction. 
An interesting observation - for which we have no concep- 
tual explanation at the moment - is that the variability 
of the magnitudes of the Fourier coefficients decreases as 
the order of the critical point increases. 

Figures [s] and [9] illustrate the CLP method (The 



leo- 



III. 1[ ) in practice. In the top part of each figure 



we have plotted on a log-log scale the norms in the left- 
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FIG. 7: Plot of logiQ |fc (3?x)fc| vs. logj^Q jfcj for the maps 
/d,(2°°>,i+3i with d = 1, 5, 20. 
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FIG. 9: Top: Plot of log^p iV^(T) vs. log^o r for the map ^ 
with G — (5°°). Bottom: Plot of the first differences between 
the vertical coordinates of adjacent points in the top figure. 



11 = 1 
n=2 



11 = 4 



11 = 5 



11 = 6 



FIG. 8: Top: Plot of log^g iV,,(r) vs. log^, r for the map Q 
with a — (1°°). Bottom: Plot of the first differences between 
the vertical coordinates of adjacent points in the top figure. 



hand side of (11) 



d_ 
dr 



(12) 



L~(T) 



as functions of r, for ry = 1, 2, . . . , 6 for the map ([3| with 
for rotation numbers (1°°) and (5°°), respectively. For 
each value of 77, the "line" consists of 400 points corre- 
sponding to 400 different values of r for which we have 
computed the corresponding norm. 

The bottom parts of Figures |8] and [9] show the behavior 
of the first differences between the vertical coordinates of 
adjacent points from the top parts of the figures. Clearly, 
the points in the top parts do not lie on exact straight 
lines but have small periodic (as functions of log r) dis- 
placements. To make this more clear in the bottom part 
of Figure |9] we have plotted the points for 77 = 6 with 
larger circles, and have connected them with lines. 

We would like to point out that the bound (111 is only 
an upper bound. However, the discrete scaling invariance 
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. T ^, , , TT , , , r , late regularly. For this range of login t, we found the 

TABLE I: Global Holder regularity k lor maps with rotation ^ r i r at / ^ /TTTn, i • 

, /,oo\ t n A r\ 1 -i- 1 • i r rate 01 clecav 01 the uomis A'„ ( T ) 12 by measuring the 

numbers a = (k } (k = 1,2,3,4,5) and critical point of or- , ^ ... ' ^—r'^ fi i •? 

slopes, K — T], irom which we computed the regularity k. 

The accuracy of these values is difficult to estimate, 

but a conservative estimate on the relative error of the 

data in Table |l] is about 3%. 

We have also computed the regularity k and the scaling 

exponent a of several maps with rotation numbers with 

the same tails of the continued fraction expansion but 

with different heads. 

Within the accuracy of our computations, the results 

did not depend on the head, which is consistent with the 

predictions of the renormalization group picture. 



der d. 


d 


(1°°) 


(2°°) 


(3°°) 


(4°°) 




1 


0.621 


0.617 


0.607 


0.596 


0.578 
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0.432 


0.427 


0.417 


0.404 


0.388 


3 


0.328 


0.324 


0.313 


0.300 


0.291 


4 


0.263 


0.260 


0.252 


0.244 


0.232 


5 


0.220 


0.217 


0.210 


0.203 


0.193 


6 


0.189 


0.186 


0.180 


0.174 


0.163 


10 


0.121 


0.120 


0.115 


0.111 


0.105 


15 


0.084 


0.082 


0.079 


0.077 


0.074 


20 


0.064 


0.063 


0.061 


0.058 


0.055 



of the Siegel disk at small scales - manifested also by the 
existence of the "periodic windows" in Figure [6]- implies 
that the leading behavior of logj^Q Njj{t) as a function of 
lo 

Sio is indeed linear, and that pi) is close to being 
saturated. 

The discrete scaling invariance at small scales allows 
for log;^o ^ri{T) as a function of log^^Q r to have small pe- 
riodic corrections superimposed on the leading linear be- 
havior. Our calculations are precise enough that these 
corrections to the leading behavior are clearly visible. 
The bottom parts of Figures |8] and [9] show exactly these 
corrections. 

The form of these small periodic corrections depends 
in a complicated way on the the behavior of the Fourier 
coefficients in the "periodic windows" in the Fourier spec- 
trum. The presence of this periodic corrections is a good 
indicator of the ranges of t which are large enough that 
the asymptotic behaviour has started to take hold, but 
small enough so that they are not dominated by the 
round-off and truncation error. In our previous works 
[131 m] , we have also found periodic corrections to the 
scaling in other conjugacies related to the regularity of 
conjugacies of other critical objects. 



D. Numerical values of the Holder regularity 

In Table |T] we give the computed values of the global 
Holder regularities of the real and imaginary parts of the 
dynamically natural parameterizations x ([s]), ([t]) of the 
boundaries of the Siegel disks. We studied maps of the 
form Q, with different values of /3, and with different 
orders d of the critical point in dD; some runs with the 
map (|3| (for which the critical point is simple) were also 
performed. 

To obtain each value in the table, we performed the 
procedure outlined in Section [ill B | for at least two maps 
of the form Q . For each map we plotted the points from 
the CLP analysis for 77 = 1, 2, . . . , 6 as in the top parts 
of Figures [8| and [9] looked at the differences between the 
vertical coordinates of adjacent points (i.e., at graphs like 
in the bottom parts of Figures [s] and |9| , and selected a 
range of values of logj^Q r for which the differences oscil- 



E. Importance of the phases of Fourier coefficients 

In Figure [6| we saw that the modulus of the Fourier 
coefficients of 3?% and Q'x decreases very approximately 



< 



C 



< 



c 



(13) 



and that these bounds come close to be saturated in- 
finitely often. This seems to be true independent of what 
the rotation number is. 

This rate of decay of the Fourier coefficients would be 



implied by 3?% and 5^ being C^, but from (13 1 we can- 



not conclude that 3?% and 5x ^'^^ even continuous (note, 
for example, that the function f{x) = X^feLi \ is 
discontinuous at x = 0). It is well-known from harmonic 
analysis that the phases of the Fourier coefficients play a 
very important role, and changing the phases of Fourier 
coefficients changes the regularity of the functions (see, 

e.g., m) 



00 ^ 



1 



k=2 



fcl/4 



|a;|3/4 



as \x\ 



while 



k=2 ' ' 

In Figure |T0| we depict the phases of the Fourier co- 
efficients of for the map /o,(i°°).i whose only critical 
point, c = 1, is simple (i.e., of order d = 1). We see 
that, for small fc, the phases have a repeated pattern. If 
we consider k G Ij := [a~^ , a~^~^] (where a = (1°°) is 
the golden mean), we see that the phase restricted to Ij 
has a pattern very similar to Ij+i, except that the lat- 
ter is reversed and amplified. Of course, since the phase 
only takes values between — tt and tt, the amplification of 
the patterns can only be carried out a finite number of 
times until the absolute values of the phases reach tt, af- 
ter which they will start "wrapping around" the interval 
[— TT, tt]. Unfortunately, to see this effect numerically, we 
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FIG. 10: Plot of the phases of the Fourier coefficients 
(vs. logj^o 1^1) for the map /o,(i^),i. 



would need hundreds of millions of Fourier coefficients, 
which at the moment is out of reach. 

Given the above observation, it is natural to study the 
distribution of the phases of the Fourier coefficients in 



an interval of self-similarity. In Figure 1 1 we present the 
histogram of the phases in the interval k e I24 U I25 = 



FIG. 12: Kolmogorov-Smirnov test 
phases in I24 U I25 ■ 



for normality of the 



a KS-test and the command qqline displays the result 
of a fitted Gaussian). The Kolmogorov-Smirnov test re- 
veals that, as expected (since the variable is an angle), 
the distribution of the phases has discrepancies with a 
Gaussian near the edges, — tt and tt. Nevertheless, there 
is a remarkably good fit away from these edges. For the 
intervals we chose, most of the data points are indeed out 
of the edges. 



F. Data on the scaling exponents 

In Table |TT] we give the values of the modulus of the 
scaling exponent a for maps of the form ^ with orders 
d = 1,..., 6, 10, 15, 20, 40, 60, 80, 100 of the critical point 
and rotation numbers cr — {k°°) with k = 1, ... ,5. We 
believe that the numerical error in these values does not 
exceed 2 in the last digit. 

In [5 , the author computed the scaling exponents for 
maps with rotation number (1°°) and critical point of 
different orders d, and suggested that the behavior of a 
for large d is 



FIG. 11: Histogram of some of the phases shown in Figure [To] 
(from J24 U J25). 



as d 



(14) 



^-26] j,^^ g^^Q^^ 170,000 phases). We note that 
the histogram is very similar to a Gaussian. This vi- 
sual impression is confirmed by using the Kolmogorov- 
Smirnov test, shown in Figure [l2j Recall that the 
Kolmogorov-Smirnov test consists in plotting the empiri- 
cal distribution versus the theoretical one (for details see, 
e.g., Chapter 7]). If indeed the empirical distribution 
was a sample of the theoretical distribution, we would 
get a set of points close the diagonal. The Kolmogorov- 
Smirnov test is available in many statistical packages (we 
used the package R, in which the command qqnorm gives 



We studied the same problem for other rotation num- 
bers and, taking advantage of the extended precision, we 
carried out the computation for rather high degrees of 
criticality (w 300), see Table |l!) In Fi gurellSl we plotted 
1/(1 — \a\) versus d for five rotation numbers. Our data 
that for high values of d the modulus of a tends to 1 for 
any rotation number. The values of the constants Aa- in 
(14 1 for rotation numbers (fc°°) with k — 1,2,3,4,5 are 
approximately 0.646, 1.168, 1.531, 1.960, 1.925, respec- 
tively (the linear regression was based on the values for 
d = 40,...,300). 
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TABLE 11: Values of \a\ for maps of the form fd,a,0 Q with different rotation numbers a and orders d of the critical point. 



d 


\^ 1 


1 


\^ 1 




(5°°) 


I 


n 741 7fl 


581 1 1 30545 


484541 091 


49463945Q 


38576Q994 

xj.ooo 1 xjiy ^lyrt 


2 


n SI 91 ^S1 074-0 


68601 3947 


607981 933 


55899367 


5968809 


Q 
O 








fi441 41 Q 

U . U 4:4: _L 4: _L C7 


fil 7QQfi4 
u.u J- / yyu^: 


4 






73801 5 


7097580 




5 


SQS71 SI 


894557 


775859 


745031 9 


795495 


Q 


Q1 9*^40 


84731 4 


804946 


7768796 


759458 


10 


Q4S0S7 


89969 


87096 


851 41 6 

XJ .OO A-^±XJ 


839375 




u.you^oo 


n Q9Q77 


u.yuool 


u.oyozoo 


U.ooDOo 


20 


0.969717 


0.94601 


0.92971 


0.919149 


0.91236 


40 


0.984450 


0.97211 


0.96356 


0.9579 


0.9544 


60 


0.989504 


0.981138 


0.975314 


0.97151 


0.96906 


80 


0.9920788 


0.98575 


0.981334 


0.97845 


0.97658 


100 


0.993638 


0.98854 


0.98499 


0.98267 


0.98116 


200 


0.996793 


0.99422 


0.99242 


0.991241 


0.99048 


300 


0.997857 


0.99613 


0.99493 


0.994140 


0.99363 




FIG. 13: Plot of 1/(1 - |a|) vs. d for rotation numbers (1°°) 
(circles), (2°°) (squares), (3°°) (diamonds), (4°°) (triangles), 
(5°°) (pluses). 



IV. CALCULATION OF THE SIEGEL RADIUS 
AND THE AREA OF THE SIEGEL DISK 

As a byproduct of our calculations we can obtain rather 
precise values of two quantities of mathematical interest: 
the area of the Siegel disk and the Siegel radius. 



A. Calculation of Siegel radius 

We note that the parameterization, x, of the boundary 
dV is related to the conjugacy /i ([2| by ([s]). Then, the 
Fourier coefficients of % satisfy |x-fcT~ l^fer s ^"-"^ A: e N, 
where hk are the Taylor coefficients of h (recall that /iq = 
and hi = 1). As shown in ^HIIT], one can get the all 
the coefficients hk by equating terms of like powers in (|2| , 



and this gives infinitely many different ways to compute 
rs. In particular, since hi — 1, we have rg — |x-i|. Since 
we also have h2 — /2/[a(a — 1)] and \a\ — 1, we obtain 
''i = |a - IIIX-2I/I/2I (where /i = a, /2, ... are the 
Taylor coefficients of the function /). Similar formulas 
for higher order terms are also available. 



B. Calculation of the area of the Siegel disk 

Since h{ra^z) is a univalent mapping from the unit disk 
to the Siegel disk, we can use the area formula [15] 



Area — 



^ k\hk\^ ^ kr 



-2k\ 



fe=l 



IX-fcl 



(15) 



For polynomials, the Siegel disk is bounded so that 
the sum in (|15[ ) is finite. This is compatible with the 
observation ( |13[ ), but it shows that the bound cannot be 
saturated very often. 

In Table we give the values of the areas of the Siegel 
disks of the map fd,a,t} Q (we believe that the error 
does not exceed 2 in the last digit). Since the series (151 
converges slowly, we computed the partial sums of the 
first Qn terms in (15 1, where Qn are the denominators 
of the best rational approximants, Pm/Qm = {k"^), to 
the rotation number a = {k°°) (cf. (|9|), and then per- 
formed Aitken extrapolation on these values. Because of 
the repeating "periodic windows" in the Fourier spectra 
(shown in Figures |4)[5)[6]|7]), these partial sums tend to 
the area of the Siegel disk geometrically, and Aitken ex- 
trapolation gives good results. In our computations we 
used 2^^ « 8 X 10^ Fourier coefficients of x- 

Clearly, the area of a Siegel disk depends on the partic- 
ular choice of the map /, i.e., is non-universal. Perhaps 
the only universal characteristic that can be extracted is 
the rate of convergence in the Aitken extrapolation, but 
we have not studied this problem in detail. 
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TABLE III: Areas of the Siegel disks for nraps of the form /d,{fe=°),i+3i (see Q) with different rotation numbers a and orders 
d of the critical point. 



d 


(1°°) 


(2°°) 


(3°°) 


(4°°) 


(5°°) 


1 


1.3603361 


1.3586530 


1.3611085 


1.3652030 


1.3693337 


9 

£i 


u.oyoDOy 




u.oyoDUO 




u.oyoD ( 


3 


0.65986 


0.65766 


0.65664 


0.6553 


0.65308 


4 


0.5190 


0.5170 


0.51550 


0.5133 


0.5104 


5 


0.4262 


0.4244 


0.4226 


0.4200 


0.4170 


6 


0.3607 


0.3591 


0.3573 


0.354 


0.3515 


10 


0.2214 


0.2203 


0.2187 


0.2163 


0.2137 


15 


0.148 


0.147 


0.146 


0.144 


0.1421 


20 


0.111 


0.110 


0.109 


0.107 


0.106 



TABLE IV: Upper bounds on the Holder regularity K^ax 
for maps of the form fd,a;p Q with rotation numbers a 
[k = 1, 2, 3, 4, 5) and critical point of order d. 



d 


(1-) 


(2°°) 


(3°°) 


(4°°) 


(5-) 


1 


0.6203 


0.6159 


0.6064 


0.5933 


0.5783 


2 


0.4324 


0.4276 


0.4175 


0.4039 


0.3890 


3 


0.3293 


0.3252 


0.3164 


0.3047 


0.2922 


4 


0.2653 


0.2618 


0.2543 


0.2443 


0.2338 


5 


0.2219 


0.2189 


0.2124 


0.2039 


0.1949 


6 


0.1906 


0.1880 


0.1823 


0.1749 


0.1670 


10 


0.1218 


0.1200 


0.1163 


0.1114 


0.1063 


15 


0.0838 


0.0826 


0.0800 


0.0766 


0.0731 


20 


0.0639 


0.0630 


0.0610 


0.0584 


0.0557 



V. AN UPPER BOUND OF THE REGULARITY 
OF THE CONJUGACY 

As pointed out in [T^l Section 8.2], one can find upper 
bounds for the regularity in terms of the scaling expo- 
nents. 

Recall that x • T — > dV conjugates (the restriction 
to dV of) the map / : dV dV to the rigid rotation 
: T ^ T : t 1-^ t + a (where a is the rotation number 
of /), namely, X ° ''a- = / ° X- Let us consider only rota- 
tion numbers of the form a — {k°°), and let the natural 
numbers Qm and the scaling exponent a be defined by 
([9]) and (10 1. Then closest returns of the iterates of € T 
of To- and the iterates of c G dV to the starting points 
and c, respectively, are governed by the scaling relations 



|/^"(c)-c| 



Cia"-|-o(a™) , 
C2\ar + o{\an 



So that we obtain that h{Cia"^) ~ C2\a\'' 
This is impossible if h is C" with 

K > Kmax = log |a|/log(J 



(16) 



and the right-hand side of ( 16 1 is not an integer. 



In Table |IV] we give the values of the upper bound to 
the regularity Kmax ( 16 1 for maps with rotation numbers 
and order of the critical point. To compute these values, 
we used the values for \a\ from Table |ll] (and the exact 



values for the rotation numbers) . We have kept only four 
digits of accuracy, although the error of these numbers 
is smaller (their relative error is the same as the relative 
error of the values of \a\). 

Clearly, within the numerical error, the values of the 
regularities from Table |l] (obtained from applying the 
CLP method) are equal to the upper bounds on the reg- 
ularity from Table IV (obtained from the scaling expo- 
nents) . 

This is in contrast with the results in [13] , where similar 
bounds based on scaling were found to be saturated by 
some conjugacies but not by the inverse conjugacy. 

For maps with highly critical points, if (16 1 holds, 
then the asymptotic behavior of |a| (14 1 implies that 
the asymptotic behavior of the limit on the regularity 
becomes 



log 1 



log cr 



|log|cr| 



as d 



(17) 



VI. CONCLUSIONS 



We have considered Siegel disks of polynomials with 
some quadratic fields and with different degrees of critical 
points. 

We have made extended precision calculations of scal- 
ing exponents and a parameterization of the boundary. 
This allows us to compute the regularity of the boundary 
using methods of harmonic analysis. 

The regularity of the boundary seems to be univer- 
sal, depend only on the tail of the continued fraction ex- 
pansion and saturate some easy bounds in terms of the 
continued fraction expansion. 

We have identified several regularities of the Fourier se- 
ries of the conjugacy. Namely, it seems that, irrespective 
of the rotation number, we have < limsup|fcx/c| < 00. 
There seems to be a regular distribution of the phases of 
the Fourier coefficients which follows a Gaussian law. 

We have also extended the results of |3] on the depen- 
dence of scaling exponents on the degree of the critical 
point to higher degrees and to other rotation numbers. 
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